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Using a covariant coalescence model, we study hadron production in relativistic heavy ion collisions 
from both soft partons in the quark-gluon plasma and hard partons in minijets. Including transverse 
flow of soft partons and independent fragmentation of minijet partons, the model is able to describe 
available experimental data on pion, kaon, and antiproton spectra. The resulting antiproton to 
pion ratio is seen to increase at low transverse momenta and reaches a value of about one at 
intermediate transverse momenta, as observed in experimental data at RHIC. A similar dependence 
of the antikaon to pion ratio on transverse momentum is obtained, but it reaches a smaller value 
at intermediate transverse momenta. At high transverse momenta, the model predicts that both 
the antiproton to pion and the antikaon to pion ratio decrease and approach those given by the 
perturbative QCD. Both collective flow effect and coalescence of minijet partons with partons in 
the quark-gluon plasma affect significantly the spectra of hadrons with intermediate transverse 
momenta. Elliptic flows of protons, Lambdas, and Omegas have also been evaluated from partons 
with elliptic flows extracted from fitting measured pion and kaon elliptic flows, and they are found 
to be consistent with available experimental data. 
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I. INTRODUCTION 



Recently, there is a renewed interest in using the par- 
ton coalescence or recombination model to understand 
hadron production from the quark-gluonplasma formed 
in relativistic heavy ion collisions [lj, |2|, |3j, U, |5|, |6| ■ Em- 
phases in these studies are, however, different from earlies 
ones based on the coalescence model such as the ALCOR 
and the MICOR || model. Instead of addressing par- 
ticle yields and their ratios, these new studies were more 
concerned with observables related to collective dynam- 
ics and production of hadrons with relatively large trans- 
verse momentum. In Ref. p] , the parton coalescence was 
used to convert the quark matter, that is formed from 
melted soft strings produced in initial soft collisions, to 
hadrons. Including parton coalescence in a multiphase 
transport model (AMPT) it was found that partonic 
effects are important for describing measured large el- 
liptic flows and narrow two-pion correlation functions at 
RHIC. In Ref. |2|, it was shown that hadron production 
based on parton coalescence is able to account for the 
qualitative difference between the observed elliptic flows 
of mesons and baryons. Based on parton recombination, 
a parton transverse momentum distribution is obtained 
in Ref. 3] from the measured pion spectrum and is then 
used to predict the kaon and proton transverse momen- 
tum spectra. The parton coalescence model is further 
found in Refs. 0, Q to be able to explain the observed 
enhancement of intermediate transverse momentum pro- 
tons and antiproton. While coalescence of partons from 
a quark-gluon plasma with high effective temperature is 
considered in Ref. 3, their coalescence with minijet par- 
tons is introduced in Ref. la as a new mechanism for 



hadronization of minijet partons. 

Besides independent fragmentations to hadrons as usu- 
ally considered, minijet partons produced in relativistic 
heavy ion collisions |l(j are allowed in Ref. t 5] to coa- 
lesce with partons from the quark-gluon plasma formed 
in relativistic heavy ion collisions, as suggested in Ref. 
6] for studying the flavor ordering of the elliptic flows 
of hadrons with intermediate transverse momenta. Since 
minijet partons have a power-law transverse momentum 
spectrum while partons in the quark-gluon plasma have 
an exponential thermal spectrum, this mechanism leads 
to an enhanced production of hadrons with intermedi- 
ate transverse momentum. Because of stronger enhance- 
ment for baryons and antibaryons than for pions from 
this hadronization mechanism, a large antiproton to pion 
ratio of about one is obtained at intermediate trans- 
verse momenta as seen in the experimental data from 
the PHENIX collaboration ^lj. However, in this study 
hadrons from hadronization of the quark-gluon plasma 
are taken to also have exponential thermal spectra ex- 
tending to all transverse momenta. Furthermore, in or- 
der to obtain a semi-analytical expression for the coales- 
cence formula, only comoving partons at zero rapidity are 
considered in evaluating the coalescence probability for 
hadron production. Although relativistic kinematics was 
used in this study, the model is not fully covariant. In the 
present paper, we relax these simplifications by using a 
covariant coalescence model and treating more generally 
parton coalescence via a Monte-Carlo method. Also, we 
include coalescence among partons from the quark-gluon 
plasma so that hadrons with low momenta are treated on 
the same footing as those with intermediate momenta. 
We show that this improved coalescence model for par- 
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tons from both the quark-gluon plasma and minijets is 
able to reproduce the experimental transverse momen- 
tum spectra of pions, antikaons, and antiprotons mea- 
sured at RHIC as well as the antiproton to pion ratio as 
functions of transverse momentum. We further predict 
the dependence of the K~ /7r~ ratio on transverse mo- 
mentum. Fitting quark elliptic flows to measured pion 
and kaon elliptic flows, the predicted elliptic flows of pro- 
tons and A are found to agree with available experimental 
data. The f2 elliptic flow is also predicted and is smaller 
than that of A as a result of smaller elliptic flow for light 
quarks than strange quarks. 

The paper is organized as follows. In Section [HI we 
describe the general formalism of a covariant coalescence 
model for mesons and baryons. How minijet and quark- 
gluon plasma partons are determined for Au+Au colli- 
sions at 200 AGeV are described in Section lull In Sec- 
tion [W] the Monte-Carlo method used for treating coa- 
lescence among partons whose numbers vary by many or- 
ders of magnitude is presented. Results for the transverse 
momentum spectra of pions, antiprotons, and antikaons 
obtained from the coalescence model are given in Section 
Ivl Ratios of the antiproton to pion and the antikaon to 
pion transverse momentum spectra are also shown. We 
further study the effect of coalescence of minijet partons 
with those in the quark-gluon plasma and also the effect 
of collective flow on the spectra and ratios of produced 
hadrons. The elliptic flow of hadrons based on that of 
quarks is also studied. Finally, we conclude in Section 
IVIIII with a summary of present work and an outlook 
about future developments and applications of the par- 
ton coalescence model. 



II. THE COALESCENCE MODEL 

Using the covariant coalescence model of Dover et at 
[l2|. the number of mesons formed from the coalescence 
of quark and antiquarks can be written as 



N M = 9m I Pi ■ da x p 2 ■ da 2 



d 3 pi d 3 p 2 
(2tt) 3 E 1 (2tt) 3 E 2 
x f Q {x 1 ;p 1 )fq(x 2 ;p 2 )f M (x 1 ,x 2 ;p 1 ,p 2 ). (1) 

In the above, da denotes an element of a space-like hy- 
persurface; and qm is the statistical factor for forming a 
colorless meson from spin 1/2 color quark and antiquark. 
For mesons considered here, i.e., 7r, p, K, and K* , the sta- 
tistical factors are g„ = gK = 1/36 and g p = gK* = 1/12. 
The functions f q (x,p) and fq(x,p) are, respectively, co- 
variant distribution functions of quarks and antiquarks in 
the phase space, and they are normalized to their num- 
bers, i.e., 



p ■ da 



d 3 i 



{2ir) 3 E 



fq,q( X iP) 



AW 



(2) 



The function fAii^i, x 2 ;pi,p 2 ) in Eq.Q is the proba- 
bility for a quark and an antiquark to form a meson. It 



describes the dynamic process of converting a quark and 
an antiquark to a bound state meson in the presence of a 
partonic matter. It depends on the overlap of the quark 
and antiquark wave functions with the wave function of 
the meson as well as the interactions of emitted virtual 
partons, which are needed for balancing the energy and 
momentum, with the partonic matter. Neglecting the 
off-shell effects and taking the wave functions of quark 
and antiquark to be plane waves, the coalescence proba- 
bility function is then simply the covariant meson Wigner 
distribution function. In Ref. |l2j| . it is parametrized by 
Gaussian functions in x\ — x 2 and p\ —p 2 . Here, we take 
it to have a uniform distribution as in Ref. jj, i.e., 



fM(xi,x 2 ;pi,p 2 ) 



9tt 



r© (A^ - (xi - x 2 ) 2 ) 



2(A X A P )3 

x6 (A 2 p - ( Pl - p 2 f + (mi - to 2 ) 2 ) , (3) 

where and A p are the covariant spatial and mo- 
mentum coalescence radii. The factors in front of theta 
functions are introduced to obtain the correct normaliza- 
tion for the meson Wigner function in the nonrelativistic 
limit. Here we use h = c = 1. 

For ultrarelativistic heavy ion collisions at RHIC, it is 
convenient to introduce rapidities variables y and rj in 
the momentum and the coordinate space, respectively, 
and they are defined by 



y = — in — 

y 2 E 



Pz 



1. t + z 
7? = — in . 

' 2 t-z 



(4) 



Using these variables, the spatial coordinate becomes 
x = (rcoshm — rsinhm — ]"t), with r = \/t 2 — z 2 
and i"t denoting, respectively, proper time and trans- 
verse coordinates; while the momentum is p = 

(tot coshy, —tot sinhy, — pt) with tot = \ nr- 



be- 
ing the transverse mass in terms of quark mass m q and 
its transverse momentum pt- The momentum volume 
element is then given by 



d3 P ,1 J2 
— = dyd p T , 

while the spatial volume element becomes 
p ■ da = ttot cosh(y — rfjdrjd 2 ^, 



(5) 



(0) 



if we adopt a hypersurface of constant longitudinal 
proper times. 

The quark and antiquark phase space distribution 
functions arc then given by 



fq,q( X 'P) 



(27T) 3 



dN„ 



titit cosh(y — rj) dr/d 2 rTdyd 2 pT ' 



(7) 



where R± is the transverse radius of the system. 

The yield of mesons from coalescence of quarks and 
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antiquarks is then given by 



N M = 9m \ dT] 1 d 2 r 1T dT] 2 d 2 r2Tdyid 2 p 1T dy 2 d 2 p2T 



dN a 



dNg 



drjid 2 TiTdyid 2 piT d^d 2 v 2 Tdy 2d 2 P2T 

X fM{xi,X-z)Px,P2). (8) 

The transverse momentum spectrum of mesons can be 
obtained from Eq. 1|12|) by multiplying the right hand side 
with 



1 = J d 2 p T 5 ( - 2) (p T - Pit - Pst), 



(9) 



and then differentiating both sides of the equation with 
respect to px- The resulting meson transverse momen- 
tum spectrum from quark and antiquark coalescence is 
given by 



= g M I dr] 1 d 2 r 1T dri 2 d 2 r 2T dy 1 d 2 p 1 T 



dN n 



dNg 



dN M 
d 2 pr 

dy 2 d P 2T ^ r?i ^2 riT ^y iC ;2p iT dr\2d 2 ?2Tdy2d 2 P2T 
xf M (x 1 ,X2]Pi,p 2 )S ( - 2 \p T - Pit - P2t)- (10) 

For quarks and antiquarks produced at central rapidi- 
ties in relativistic heavy ion collisions, it is reasonable to 
assume that their longitudinal momentum distributions 
are boost-invariant, i.e., independent of rapidity. Fur- 
thermore, they are expected to satisfy the Bjorken cor- 
relation of equal spatial and momentum rapidities, i.e., 
rj = y. The quark and antiquark phase space distribution 
functions in the rapidity range Ay can then be expressed 
as 



dN„ 



5(r)-y) dN q 



drjd 2 rTdyd 2 pT 



Ay G? 2 rT<i 2 PT 



i2/l<Ay/2 



(11) 



This leads to the following meson transverse momentum 
spectrum from coalescence of quarks and antiquarks 



dN 



M 



d 2 PT 



9M 

(Ay) 2 



dN n 



d 2 r 1T d 2 r 2T d 2 p 1T d 2 p 2T 



dN, 



e? 2 riTG? 2 PiT 



ll/2|<A a /2 



\yi\<Ay/2d 2 Y 2T d 2 p 2T 

driidy d r]2dy2S(rj 1 - yi)8(r]2 - y 2 ) 
x fM(xi 1 x 2 ;p 1 ,p 2 )5 {2) (v T - Pit - P 2 t)- (12) 

The above result can be simplified if there is no corre- 
lation between parton transverse momenta and positions, 
such as in the absence of collective transverse flow, and if 
partons are also uniformly distributed in the transverse 
space. In this case, the quark and antiquark distributions 
only depend on transverse momentum, i.e., 



dN n 



cPttgPpt 



1 dN q<q 



\y\<by/2 TlR 2 , d 2 p T 



\y\<&y/z 



(13) 



Considering small rapidity range such as Ay < 1, we 
have 



(xi-x 2 ) 2 = 2r 2 [l - cosh(ryi - 772)] - (r iT - r 2T ) 



-r 2 {m - rj 2 ) 2 
-(ri- r 2 ) 2 , 



(riT - r 2T ) 



(14) 



and 



(Pi ~ P2) = m 1T + m 2T - 2m 1T m 2 T cosh(yi - y 2 ) 

^2 



(PIT — P2T) 

(miT - m 2T ) 2 - (Pit - P2t) 



(15) 



After carrying out the integrals in transverse space, the 
meson transverse momentum spectra is then give by 



dN 



Gtt 



d 2 pi 



9m 



dN n 



d?p 



IT 



dN- n 



\yi\<&y/2d 2 p2T 



d 2 Pnd 2 p 2 T 



y 2 \<Ay/2 



>.5 {2 \p T - PlT - p 2T )e (a 2 - (p 1T - p 2T ) 2 

-[(miT - m 2T ) 2 - (mi - m 2 ) 2 ]) . 



(16) 



This result for meson transverse momentum spectra re- 
duces to that used in our previous schematic study if 
we assume that only comoving quarks and antiquarks can 
coalesce to hadrons and replace the arguments of theta 
functions in Ea. H16|) by A 2 minus quark and antiquark 
relative momentum in meson rest frame. 

Since we will include collective transverse flow of par- 
tons in the quark-gluon plasma, Ea. l|12JI will be used in 
the following study. To generalize the results for mesons 
to formation of baryons and antibaryons from the parton 
distribution functions, we take the baryon coalescence 
probability function as 

F B (xi,X2,x 3 ;p 1 ,p2,p 3 ) 

2 1 

9tt „ / . , 1 



x9 [A 2 --( Pl - P 2) 



xe^A 2 -i[( Pl+J32 -2p 3 ) 2 
-(m 1 +m 2 -2m 3 ) 2 ]), (17) 

where we have taken for simplicity the same space and 
momentum coalescence radii for the relative Jacobi coor- 
dinates among three quarks. 

For boost-invariant dynamics with Bjorken spatial and 
momentum rapidities correlation, we obtain following 
baryon transverse momentum spectrum from coalescence 
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of three quarks: 
dN B cjb 



r 3 

J Y[d 2 r iT d 2 piT- t 



dN n 



d 2 r iT d 2 p iT \ yi \<^ y /2 

i—l 



(18) 



In the above, qb is the statistical factor for formation of a 
baryon from three quarks. For baryons and antibaryons 
considered in present study, i.e., p, A, p and A, the sta- 
tistical factors are g p — g p — 1/108 and g& = g& = 1/54. 
The above formula can also be used for antibaryons by 
replacing quark momentum spectra by the momentum 
spectra of antiquarks. 



III. PARTON DISTRIBUTIONS 

We consider central Au+Au collisions (0 — 10%) at 
200 AGeV, available at RHIC. In these collisions, hard 
processes between initial nucleons lead to production of 
minijets with large transverse momentum. Also, a quark- 
gluon plasma is expected to be formed from soft processes 
in the collisions. In this section, both parton momentum 
spectra in minijets and in the quark-gluon plasma are 
introduced. 



A. minijet partons 

Partons at high transverse momenta (usually greater 
than 2 GeV) are mainly from the minijets produced 
in initial hard collisions among nucleons. The trans- 
verse momentum distribution of minijet partons in the 
midrapidity can be obtained from an improved pertur- 
bative QCD calculation It is given by dNj et /d 2 pT = 
1 /<Jtot W do jet /d 2 pr in terms of cr^~ t 10 corresponding to 
the total cross section at central 10% of the collisions 
and the jet production cross section from nucleus-nucleus 
collisions, 



da-. 



jet 



J d 2 b d 2 rt Au (r)t Au (b-r) 
x ^ Jdx a dx b d 2 k aT d 2 k bT g(k aT )g(k hT ) 

ab 



X fa/Au(x a ,Q )fb/Au(x b ,Q ) 

S r „ £ „,d<r ab 
x — d(s + t + u) — — . 

7T dt 



(19) 



In the above, t A (r) is the thickness function of Au at 
transverse radius r given by integrating the nuclear den- 
sity distribution along the longitudinal direction. The 



parton distribution function in a nucleon in the nu- 
cleus Au is denoted by f a /Au( x >Q 2 ) including transverse 
smearing g(ki-). The cross section da ab /dt is the parton 
scattering cross section. Kinematic details and a system- 
atic analysis of pp collisions can be found in Ref. [lflj . 
Using the GRV94 LO result for the PDFpJ and the 
KKP fragmentation function from Ref. |15|. measured 
data in the reaction pp — > ir°X at yfs = 200 GeV can be 
reproduced with Q = 0.75p T and = 2 GeV 2 . 

In heavy ion collisions at RHIC, minijet partons are 
expected to lose energy by radiating soft partons as they 
traverse through the quark-gluon plasma. This effect can 
be taken into account by lowering their transverse mo- 
menta by the energy loss A_E, which depends on both 
the parton energy E and an effective opacity parameter 
L/X according to the GLV model [la |. 




3 4 5 

P T (GeV) 

FIG. 1: Transverse momentum distributions of partons at 
hadronization in Au+Au collisions at i/s = 200 AGeV for 
gluons (short-dashed curve), u and d (solid curve), u and d 
(dashed curve) as well as s and s (dash-dotted curve) quarks. 
Minijet partons have transverse momenta greater than 2 GeV, 
while partons from the quark-gluon plasma have transverse 
momenta below 2 GeV. 

Taking as momentum cutoff po = 2 GeV/c for mini- 
jet partons and using an effective opacity L/X = 3.5 as 
extracted from fitting the spectrum of high transverse 
momentum pions measured at RHIC ^a, , the trans- 
verse momentum spectra of minijet partons at midrapid- 
ity (y = 0) in central Au+Au collisions at yfs = 200 
AGeV are shown in Fig. ^ f° r gluons (short-dashed 
curve), u and d (solid curve), u and d (dashed curve) 
as well as s and s (dash-dotted curve) quarks. These 
spectra can be parametrized as 



rfVjc 

d 2 Pi 



A 



D 



B+p T 



(20) 



Values for the parameters A, B, and n for gluons, light 
and strange quarks and anitquarks are given in Table [I] 
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For later calculation of parton coalescence probability, 
masses of minijet partons are taken to be the current 
quark masses, i.e., 10 MeV for light quarks and 175 MeV 
for strange quarks. Also, we assume that the rapidity 
distribution of minijet partons is uniform in the rapidity 
range of y € (—0.5, +0.5) considered in the present study. 



TABLE I: Parameters for minijet parton distributions given 
in Ea. lffiljl at midrapidity from Au+Au at tfs = 200 GeV. 





A[l/GeV 2 ] 


B[GeV] 


n 


g 


3.2xl0 4 


0.5 


7.1 


u,d 


9.8xl0 3 


0.5 


6.8 


u, d 


1.9xl0 4 


0.5 


7.5 


s,s 


6.5xl0 3 


0.5 


7.4 



B. the quark-gluon plasma 

For partons in the quark-gluon plasma, their trans- 
verse momentum spectra are taken to have an expo- 
nential form. For their longitudinal momentum distri- 
bution, we assume that they are boost-invariant, i.e., 
they have a uniform rapidity distribution in the range 
y G (—0.5, +0.5). To take into account collective flow of 
quark-gluon plasma, these partons are boosted by a flow 
velocity vt = /3o( rr r/ R±), depending on their transverse 
radial positions rr. Here, R± is the transverse size of 
the quark-gluon plasma at hadronization, and /3q is the 
collective flow velocity of the quark-gluon plasma and is 
taken to be 0.5c. In this case, the light quarks and anti- 
quarks transverse momentum spectra are given by 



d 2 rxd 2 PT 



9g.,gTm T 
(2tt) 3 



exp 



7 T (m T - p T • v T + fi q ) 



T 



,(21) 



where g q = g q = 6 are the spin-color degeneracy of light 
quarks and antiquarks, and the minus and plus signs are 
for quarks and antiquarks, respectively. The slope pa- 
rameter T is taken to be T — 170 MeV, consistent with 
the phase transition temperature (T ~ 165 — 185 MeV) 
from lattice QCD calculations. Masses of light quarks 
and antiquarks are taken to be m q = m q = 300 MeV, 
similar to the masses of constituent quarks due to pos- 
sible nonperturbative effects in the quark-gluon plasma 
near hadronization |18| . For the quark chemical potential 
fx q , we use a value of \i q = 10 MeV to give a light anti- 
quark to quark ratio of 0.89, which would then lead to an 
antiproton to proton ratio of about (0.89) 3 = 0.7, con- 
sistent with the observed ratio at midrapidity in heavy 
ion collisions at RHIC. The transverse flow effect i s taken 
into account through the factor 7t = 1/ \J\ — v\. 

The momentum spectra for strange quarks and anti- 
quarks are similar to Ea. (|21|) with m q replaced by the 



strange quark mass m s — m s = 475 MeV and fi q replaced 
by /j, q — fi s . The strange chemical potential [i s is taken to 
have the same value as p q in order to have same strange 
and antistrange quarks numbers. The resulting strange 
quark to light quark ratio is then about 0.27. Including 
contribution from decays of p and K* , this would then 
give a final K~/K+ = 0.89 and K~ /<k~ = 0.24, com- 
parable to experimental data at midrapidity from RHIC. 
Eq. (|2"T)l also applies to gluons after replacing g q by the 
gluon spin-color degeneracy g g — 16 and dropping the 
chemical potential. For the gluon mass, we take it to be 
similar to that of light quarks in order to take into ac- 
count nonperturbative effects in the quark-gluon plasma. 

Although the shape of quark-gluon plasma is nearly 
cylindrically symmetric, its parton momentum distribu- 
tion is azimuthally asymmetric as shown by the finite 
elliptic flow of hadrons observed in experiments. Since 
we are interested in hadron spectra that are integrated 
over the azimuthal angle, including such asymmetry is 
not expected to change much our results. Nonetheless, 
the Monte-Carlo method to be described in Section lTvl to 
evaluate the coalescence formula can be easily extended 
to take into account such effects in order to study the 
relation between the elliptic flow of partons and those of 
hadrons as suggested in Ref. Q and to be studied in 
Section IVTT1 

The quark-gluon plasma is further assumed to have a 
transverse radius of R± = 8.3 fm at proper time r = 4 
fm, corresponding to a volume of V = 900 fm 3 . Positions 
of partons in the transverse direction are taken to have 
a uniform distribution. Their longitudinal positions are 
then determined by z — t sinh y, as we have assumed that 
rj = y due to assumed Bjorken correlation. The result- 
ing total transverse energy per unit rapidity from both 
the expanding quark-gluon plasma and minijet partons is 
about 590 GeV and is consistent with that measured by 
the PHENIX collaboration [n|. Most of this tr ansverse 
energy comes from soft QGP partons as the contribu- 
tion of minijet partons is only about 10%. The resulting 
parton density is about /9 par ton ~ 1 fm -3 . The thermal 
parton spectra below the momentum cutoff po = 2 GeV 
from the quark-gluon plasma including the collective flow 
effect are shown in Fig. \l\ior gluons (short-dashed curve), 
u and d (solid curve), u and d (dashed curve) as well as 
s and s (dash-dotted curve) quarks. Because of scat- 
tering of minijet partons with thermal partons as they 
traverse the quark-gluon plasma, those with momentum 
around po are expected to be thermalized with partons 
in the quark-gluon plasma, leading to a smooth spectrum 
around po- In the present study, we neglect this effect. 



IV. THE MONTE-CARLO METHOD 

In our previous study [5|, only partons at midrapid- 
ity iu = 0) are considered. Furthermore, only comoving 
partons, i.e., partons with momenta in the same trans- 
verse direction are allowed to coalesce to hadrons. In 
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these limits, the coalescence formula is reduced to a one- 
dimensional integral for mesons and a two-dimensional 
integral for baryons. In the present study, we do 
not introduce these simplifications. Instead, the multi- 
dimensional integrals in the coalescence formula, given 
by Eos. (|12f) and Ijl8(l are evaluated by the Monte-Carlo 
method via test particles. Specifically, we introduce a 
large number of test partons with uniform momentum 
distribution. To take into account the large difference 
between numbers of thermal and minijet partons, a test 
parton with momentum px is given a probability that 
is proportional to the parton momentum distribution, 
e.g., dN q /d 2 pT for quarks, with the proportional con- 
stant determined by requiring that the sum of all parton 
probabilities is equal to the parton number. With test 
partons, the coalescence formulas, Eas. (|12fl and (|18fl . are 
rewritten as 



dN M 
d 2 p T 



9M 



^p ? (OA(i)^ 2) (px-PiT- PiT ) 



1,3 



x fM(xi,Xj;pi,pj). 



(22) 



and 



dN, 



d 2 pr 



B \ - 



P g (i)P q (j)P q (k) 

X 5 (2) (p T - PiT - PjT - PfcT?) 

x f B (xi,Xj,Xk;pi,Pj,Pk)- 



(23) 



In the above, P q {i) and Pq(j) are probabilities carried by 
ith test quark and jth test antiquark. 

The Monte-Carlo method introduced here allows us to 
treat the coalescence of low momentum partons on the 
same footing as that of high momentum ones. We find 
that despite a decrease of eight orders-of-magnitude in 
real particle spectra, about equal numbers of test hadrons 
are formed at all momenta. 



V. HADRON TRANSVERSE MOMENTUM 
SPECTRA 

In this section, we show results for the transverse mo- 
mentum spectra of pions, antiprotons, and antikaons us- 
ing the model described in previous sections. For the 
coalescence contribution, we first take into account the 
effects due to gluons in the quark-gluon plasma and mini- 
jets by converting them to quarks and antiquark pairs 
with probabilities according to the flavor compositions 
in the quark-gluon plasma, as assumed in the ALCOR 
model 0- For both mesons and baryons, we include not 
only coalescence of hard and soft partons as in Rcf. 5] 
but also that among soft as well as hard partons. Fur- 
thermore, we include stable hadrons such as pion, nu- 
cleon (antinucleon) , and kaon (antikaon) as well as un- 
stable resonances such as p, A, and K* . Since the co- 
alescence model can be viewed as formation of bound 
states from interacting particles with energy mismatch 



balanced by other particles in the system, neglecting such 
off-shell effects is reasonable if the energy mismatch is 
small. The model is thus applicable for rho and nu- 
cleon (antinucleon) production when we take into ac- 
count quark masses. For other hadrons, the coalescence 
probability is expected to be reduced as a result of en- 
ergy mismatch. The reduction factor can in principle 
be estimated by evaluating the transition probability in 
the presence of other particles. Since the coalescence 
radii and A p are treated as parameters in our study, 
the off-shell effects can be phenomenologically taken into 
account by using different coalescence radii for different 
hadrons. However, for simplicity we use same spatial coa- 
lescence radius A x — 0.85 fm and momentum coalescence 
radius A p = 0.24 for both mesons and baryons. We note 
that there are more p and K* than pion and kaon pro- 
duced from parton coalescence. The pion, proton, and 
kaon transverse momentum spectra from parton coales- 
cence shown below include contributions from decays of 
rho, A"*, and A. 

We have also included contributions to hadron produc- 
tion from minijet fragmentations. These are obtained us- 
ing the KKP fragmentation function , which has been 
shown to reproduce measured high transverse momentum 
particles at RHIC. Explicitly, hadron momentum spectra 
are related to minijet parton momentum spectra by 



dN 



d 2 p 



had 



=£/* 



jet 



dN Dhad/jct^Q 2 ) 



'd 2 



Pjct 



(24) 



where z = Pj c t/Pha,d is the fraction of minijet momentum 
carried by the formed hadron and Q = £>had/2z is the mo- 
mentum scale for the fragmentation process. The KKP 
fragmentation function is denoted by D had /j et (z, Q 2 ). 

Before we show the results for hadron spectra and ellip- 
tic flows, we would like to point out that the coalescence 
model as formulated here is applicable if the number of 
hadrons produced is much less than the parton numbers. 
In this respect, results for hadrons with momenta above 
about 1 GeV/c, which account for about 5% of all par- 
tons, are reliable. For lower momentum hadrons, one 
needs to impose conservation of parton number in con- 
verting them to hadrons via coalescence. This correction 
has not been included in present study. 



A. pion transverse momentum spectrum 

In Fig. [3 we show the transverse momentum spec- 
trum of pions formed directly from parton coalescence 
(dashed curve). Pions from fragmentations of minijet 
partons are shown by dash-dotted curve. It is seen that 
pions from parton coalescence dominate low transverse 
momenta while those from minijet fragmentations are im- 
portant at high transverse momenta. The two contribu- 
tions have a similar magnitude at transverse momentum 
of about 3 GeV. Also shown in the figure is the total 
pion transverse momentum spectrum from the two con- 



7 



tributions (solid curve). Compared with measured spec- 
trum from the PHENIX Collaboration 2fj] (filled circles), 
the predicted spectrum of directly produced pion at low 
transverse momenta is much below experimental data. 

10 3 p . 1 . 1 . 1 . 1 . 1 . 1 . 1 . a 




012345678 
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FIG. 2: Pion transverse momentum spectra from Au+Au col- 
lisions at y/s — 200 AGeV: direct pion production from parton 
coalescence (dashed curve); pions from minijet fragmentations 
(dash-dotted curve); and sum from the above two contribu- 
tions (solid curve). Parametrized experimental tt data [2(J 
are shown by filled circles. 

Since p, K* , and A decay to pions, we have also in- 
cluded their contributions to the pion transverse momen- 
tum spectrum. It is found that rho meson decays con- 
tribute significantly to the pion spectrum at low trans- 
verse momenta and bring the final pion spectrum, shown 
by the solid curve in Fig. |3J in good agreement with ex- 
perimental data. Also shown in this figure by the dashed 
curve is the pion spectrum without including contribu- 
tion from coalescence of minijet partons with partons 
in the quark-gluon plasma. It underestimates the pion 
spectrum at intermediate transverse momenta around 3.5 
GeV/c. The contribution from coalescence of partons 
from minijets with those from the quark-gluon plasma is 
more clearly seen from the ratio of pion spectra with and 
without this contribution, shown in the inset of Fig. |3 
which is more than a factor of 2 at transverse momenta 
around 4 GeV/c. Our results thus confirm our previous 
results in Ref.[5j based on a parametrized pion spectrum 
from hadronization of the quark-gluon plasma. 

B. antiproton transverse momentum spectrum 

In Fig. 01 we show the antiproton spectrum including 
those from decays of A for Au+Au collisions at yfs = 200 
AGeV. Results for both with (solid curve) and without 
(dashed curve) contributions from coalescence of minijet 
partons with partons from the quark-gluon plasma are 
shown. They include contributions from fragmentations 
of minijet partons (dash-dotted curve). The contribu- 
tion from coalescence of soft and hard partons becomes 
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FIG. 3: Pion transverse momentum spectra from Au+Au 
collisions at y/s = 200 AGeV with (solid curve) and with- 
out (dashed curve) contributions from coalescence of minijet 
partons with the quark-gluon plasma partons. Parametrized 
experimental 7r° data |20| are shown by filled circles. Ratio of 
the two spectra is given in the inset. 
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FIG. 4: Transverse momentum spectra of antiprotons from 
Au+Au collisions yfs — 200 AGeV. The dashed curve in- 
cludes hadrons only from coalescence of partons in the quark- 
gluon plasma and from independent fragmentations of minijet 
partons (dash-dotted curve). Adding also hadrons from co- 
alescence of minijet partons with partons in the quark-gluon 
plasma gives the solid curve. Ratio of the solid to the dashed 
curve is given in the inset. Experimental p~ data |2l| at 130 
AGeV are shown by filled circles. 



important when the transverse momentum is above 3 
GeV/c, which is somewhat higher than in the case of pion 
transverse momentum spectrum. Since there are no pub- 
lished experimental data for antiproton transverse mo- 
mentum spectrum from Au+Au collisions at 200 AGeV, 
we compare our predictions with the experimental data 
from the PHENIX collaboration for Au+Au collisions at 
y/s = 130 AGeV (2JJ , shown by filled squares for trans- 
verse momenta below 3 GeV/c. Both predictions with 
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and without contributions from coalescence of minijet 
partons with partons from the quark-gluon plasma are 
comparable to the experimental data. Shown in the in- 
set of this figure is their ratio as a function of transverse 
momentum. It is seen that the ratio is close to a factor 
of 5 at transverse momenta around 5.5 GeV/c. The con- 
tribution from coalescence of minijet partons with those 
from the quark-gluon plasma is thus more important for 
antiprotons than for pions. It would be of great interest 
to have experimental data for antiprotons at such high 
transverse momenta to verify our predictions. 



C. antiproton to pion ratio 



with s+h 

without s+ 


h 


Au+Au@200AGeV 
(central) 






_\ 

\ 


^ \ 

N \ 
\ \ 
\ \ 
\ \ 



012345678 

P T (GeV) 

FIG. 5: Antiproton to pion ratio from Au+Au collisions at 
y/s = 200 AGeV. Solid and dashed curves are, respectively, 
results with and without contributions to antiprotons from co- 
alescence of minijet partons with those from the quar k-gluon 
plasma. Filled squares are the experimental data rFp| . 

The antiproton to pion ratio is shown in Fig. [S] as a 
function of transverse momentum. The solid curve is the 
result including contributions from both parton coales- 
cence and minijet fragmentations. The ratio increases 
with transverse momentum up to about 3 GeV/c and 
decreases with further increasing transverse momentum 
as in the experimental data shown by filled squares. 
At high transverse momenta, the antiproton to pion ra- 
tio becomes very small as it is largely determined by the 
results from pQCD [22}. Neglecting the contribution to 
aniproton production from coalescence of minijet partons 
with partons from the quark-gluon plasma reduces the 
antiproton to pion ratio at transverse momenta above 
2.5 GeV as shown by the dashed curve in Fig. [5j 

The enhanced antiproton to pion ratio at intermedi- 
ate transverse momenta was previously attributed to an- 
tiproton production from the baryon junctions in inci- 
dent nucleons [23|. The possibility of enhanced baryon 
to pion ratio due to parton coalescence was suggested in 
Ref. 0]. Using a parton distribution function that is 



fitted to measured pion transverse momentum spectrum, 
a parton recombination model similar to the coalescence 
model indeed leads to a large antiproton to pion ratio 
at intermediate transverse momenta 0. In Ref. Q, the 
antiproton to pion anomaly is explained by the recombi- 
nation of partons from a quark-gluon plasma with a high 
effective temperature. Our model further introduces co- 
alescence of minijet partons with partons in the quark- 
gluon plasma. This makes it possible to account for both 
the large antiproton to pion ratio at intermediate trans- 
verse momenta and its behavior at low transverse mo- 
menta. 

In our previous work [5j, an increasing antiproton to 
pion ratio at low transverse momenta was also seen, but 
it was obtained by using different inverse slope parame- 
ters for antiproton and pion transverse momentum spec- 
tra. In contrast, results obtained in the present work are 
due to coalescence of soft partons from the quark-gluon 
plasma. 

D. kaon transverse momentum spectrum and kaon 
to pion ratio 
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FIG. 6: Transverse momentum spectra of K~ from Au+Au 
collisions at yfs = 200 AGeV. The dashed curve includes 
contributions only from coalescence of partons in the quark- 
gluon plasma and hadrons from independent fragmentations 
of minijet partons (dash-dotted curve). Adding also hadrons 
from coalescence of minijet partons with partons in the quark- 
gluon plasma gives the solid curve. Ratio of the solid to the 
dashed curve is given in the inset. Experimental K~ data |2l| 
at 130 AGeV are shown by filled circles. 

In Fig. |HJ we show the antikaon spectrum includ- 
ing those from decays of K*~ for Au+Au collisions at 
^fs = 200 AGeV. Results including antikaon production 
from coalescence of soft partons from the quark-gluon 
plasma are shown by dashed curve while those including 
also coalescence between minijet partons with the quark- 
gluon plasmas are shown by the solid curve. Contri- 
butions from minijet fragmentations (dash-dotted curve) 
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are also included. As for antiprotons, we compare these 
predictions with the experimental data shown by filled 
squares from the PHENIX collaboration for Au+Au col- 
lisions at y/s = f30 AGeV |2l|, as there are no published 
data at s/s = 200 AGeV. For the limited data below 2 
GeV/c, the coalescence model reproduces them very well 
without the contribution from coalescence of minijet par- 
tons with partons from the quark-gluon plasma as the 
latter becomes important when the transverse momen- 
tum is above 2.5 GeV/c. Shown in the inset of Fig. [f)|is 
the ratio of the predictions with and without the contri- 
bution from soft and hard parton coalescence as a func- 
tion of transverse momentum. The ratio reaches a factor 
of close to 3 at transverse momentum around 4 GeV/c. 
The contribution to antikaons from coalescence of mini- 
jet partons with those from the quark-gluon plasma is 
thus comparable to that for pions. 
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FIG. 7: Antikaon to pion ratio from Au+Au collisions at 
y/s = 200 AGeV. Solid and dashed curves are, respectively, 
results with and without contributions to antikaon production 
from coalescence of minijet partons with those from the quark- 
gluon plasma. 

The antikaon to pion ratio is shown in Fig. as a 
function of transverse momentum. The solid curve is the 
result including contributions from both parton coales- 
cence and minijet fragmentations. For transverse mo- 
menta below about 2 GeV/c, this ratio is similar to the 
antiproton to pion ratio except that its value is smaller. 
At higher transverse momenta, the antikaon to pion ra- 
tio decreases slightly and reaches the value predicted by 
pQCD at high transverse momenta [24|], which gives a 
larger antikaon to pion ratio than the antiproton to pion 
ratio. Results without contribution to antikaon produc- 
tion from coalescence of minijet partons with partons 
from the quark-gluon plasma are given by the dashed 
curve, and it gives a smaller antikaon to pion ratio at 
intermediate transverse momenta compared to that with 
this contribution. Our results thus demonstrate again 
the importance of antikaon production from coalescence 
of partons from minijets and quark-gluon plasma. 
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VI. TRANSVERSE FLOW EFFECT 
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FIG. 8: Pion (upper panel) and antiproton (lower panel) 
transverse momentum spectra from coalescence of partons 
in the quark-gluon plasma with (solid curves) and without 
(dashed curves) collective transverse flow in the quark-gluon 
plasma. Ratio of these two results are shown in the insets. 



It is interesting to study the effect of transverse flow of 
quark-gluon plasma on the transverse momentum spectra 
of produced hadrons. In fact, the increase of antiproton 
to pion ratio at low transverse momenta was believed to 
be due to a stronger transverse flow effect on protons than 
on pions [25l |. Such effect was assumed in our previous 
study by using a larger inverse slope parameter for the 
antiproton transverse momentum spectrum than that for 
pions. In the present model, soft hadrons are produced 
from coalescence of partons in the quark-gluon plasma, 
which is given a collective flow velocity of 0.5c, and the 
resulting antiproton to pion ratio at low transverse mo- 
mentum is found to increase with transverse momentum. 
To see if this is due to the collective flow effect intro- 
duced in the model, we have repeated the calculations 
without transverse flow effect on antiprotons. These re- 
sults, which include only contribution to antiproton pro- 
duction from coalescence of partons in the quark-gluon 
plasma, are shown in Fig. [SI by dashed curves for pions 
(upper panel) and antiprotons (lower panel). Compared 
with results with a collective flow velocity of 0.5c, shown 
by solid curves in Fig. [SJ collective flow affects the pion 
and antiproton spectra mainly at transverse momenta 
above 1.5 GeV/c, and its effect is stronger for antipro- 
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tons than for pions. The reason that for low transverse 
momenta the inverse slope parameter for pions remains 
smaller than that for antiprotons in the absence of col- 
lective flow is due to the fact that most low transverse 
momentum pions are from decays of rho mesons, which 
gives a smaller effective slope parameter for pions than 
that of directly produced pions. Our results thus demon- 
strate that the increase seen in the antiproton to pion 
ratio at low transverse momenta is not necessarily due to 
the collective flow of quark-gluon plasma. 
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FIG. 9: Antiproton to pion ratio with (solid curve) and with- 
out (dashed curve) collective flow in the quark-gluon plasma. 
Experimental data are shown by filled squares. 

Since collective flow affects hadron spectra at inter- 
mediate transverse momenta as shown in FiglHI and its 
insets, where the ratio of pion or antiproton spectrum 
obtained with and without collective flow in the quark- 
gluon plasma is shown, it is expected to have a large 
effect on the antiproton to pion ratio at these momenta. 
This is shown in FigEI where the antiproton to pion 
ratio is given for both with (solid curve) and without 
(dashed curve) collective flow effect on protons. It is 
seen that collective flow enhances significantly this ratio 
for transverse momenta between 2 and 5 GeV/c. The 
collective flow effect on intermediate transverse momen- 
tum antiprotons is thus as strong as the effect due to 
coalescence of partons from minijets with those from the 
quark-gluon plasma, shown in Fig. |SJ To confirm the 
mechanism for antiproton production from coalescence 
of minijet and quark-gluon plasma partons, it is thus im- 
portant to have a quantitative understanding of collective 
flow in the quark-gluon plasma. 



VII. ELLIPTIC FLOWS 

The parton coalescence model based on the test parti- 
cle Monte-Carlo method is also useful for studying other 
observables at RHIC. One problem that can be addressed 
with this model is to study how elliptic flows of hadrons 



are related to that of partons prl |27|. We carry out such 
a study following the idea of Ref. where the parton 
transverse momentum distribution is extracted from fit- 
ting the pion transverse momentum spectrum using the 
quark recombination or coalescence model and is then 
used predict that of antiprotons. Here, we extract the 
elliptic flows of light and strange quarks by fitting the 
measured pion and kaon elliptic flows based on the co- 
alescence model. The resulting quark elliptic flows are 
then used to predict the elliptic flows of protons, A, f2, 
and phi mesons. 

The elliptic flow is a measure of the anisotropy of par- 
ticle transverse momentum spectrum, i.e., 



i'2 



Px -pi 



,pI 



■pi 



(25) 



where the transverse axes x and y are, respectively, in 
and out of the reaction plane. Including transverse mo- 
mentum anisotropy, the quark transverse momentum dis- 
tribution is given by 



dN n 



dN„ 



dN„ 



[1 + v 2 cos(2(f>)} , (26) 



d 2 pT p^dp^dij) p^dpT 
where </> is the azimuthal angle in the transverse plane. 




FIG. 10: Elliptic flows of pions (dashed curve) and protons 
(solid curve) as functions of transverse momentum. Elliptic 
flow of light quarks and antiquarks is given by dotted curve. 
Experimental data 28] are shown by filled squares for pions 
and filled circles for protons. 

In Fig. 1101 we show by dotted curve the elliptic flows 
of light quarks and antiquarks together with the pion el- 
liptic flow shown by dashed curve, which is supposed to 
reproduce the measured pion elliptic flow given by filled 
squares [2^. The predicted proton elliptic flow obtained 
from the quark and antiquark elliptic flows is then given 
by the dashed curve and is seen to agree with that mea- 
sured in experiments shown by filled circles p8[ . 

The elliptic flow of strange hadrons are shown in Fig. 
1111 The strange quark and antiquark elliptic flow (dot- 
ted curve) used to fit the measured kaon elliptic flow 
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FIG. 11: Elliptic flows of kaons (solid curve), phi mesons 
(dashed curve with stars), A (dashed curve), and Q (dash- 
dotted curve). Elliptic flow of strange quarks and antiquarks 
is given by dotted curve. Experimental data |2SI are shown 
by filled squares for kaons and circles for A. 

shown by filled squares for data and dashed curve 
from the coalescence model. The predicted A elliptic flow 
shown by dashed curve is seen to agree with the avail- 
able experimental data (2^| given by filled circles. We 
have also predicted the phi meson (dashed curve with 
stars) and fl elliptic flows (dash-dotted curve), which are 
smaller than the kaon and A elliptic flow, respectively, 
as a result of smaller (about 30%) strange quark elliptic 
flows than that of light quarks. 

For all hadron considered here, decrease of their elliptic 
flows at high transverse momenta is due to our assump- 
tion that high transverse momentum minijet partons and 
hadrons from minijet fragmentations have vanishing el- 
liptic flows. Of course, if high transverse momentum 
minijet partons also develop elliptic flow as suggested in 
Ref. 30], our results on hadron elliptic flows at high 
transverse momenta should be modified. 



VIII. SUMMARY AND OUTLOOK 

In this paper, we have studied the hadronization of 
quark-gluon plasma and minijet partons produced in rel- 
ativistic heavy ion collisions in terms of the parton coa- 
lescence model. The momentum spectra of partons in 
the quark-gluon plasma is taken to have an exponen- 
tial form with inverse slope parameter similar to the 
phase transition temperature, while partons in the mini- 
jets have power-law spectra. The flavor compositions in 
the quark-gluon plasma is determined by the experimen- 
tally measured antiproton to proton ratio and strange to 
nonstrange particle ratios. The volume of quark-gluon 
plasma is then fixed by the total transverse energy mea- 
sured in experiments. A collective flow is introduced in 
the quark-gluon plasma with a flow velocity comparable 
to that extracted from experiments. To take into account 



the vast difference in the magnitude of the minijet par- 
ton transverse momentum spectrum and that of partons 
in the quark-gluon plasma, a test particle Monte-Carlo 
method has been introduced to efficiently evaluate the co- 
alescence probability of partons. Both soft partons from 
the quark-gluon plasma and hard partons from minijets 
are used in the coalescence model to produce hadrons 
such as pions, kaons (antikaons), rho mesons, if*'s, nu- 
cleons (antinucleons), and A (A) resonances. Specifi- 
cally, we have included coalescence of hard minijet par- 
tons with soft partons besides that among soft and hard 
partons. 

The resulting pion, antikaon, and antiproton spectra 
are seen to agree with available experimental data from 
RHIC. For pions, contributions from rho decays are im- 
portant in explaining the measured transverse momen- 
tum spectrum below 2 GeV/c. For intermediate trans- 
verse momentum spectra between 2 and 5 GeV/c, in- 
cluding contributions from coalescence of minijet partons 
with those from the quark-gluon plasma are important, 
leading to a factor of more than two enhancement com- 
pared to results without this contribution. We have also 
compared the transverse momentum dependence of an- 
tiproton to pion ratio to the experimental data. Results 
from the coalescence model are found to agree quite well 
with the available data, i.e., it increases with transverse 
momentum and reaches a value of about one at transverse 
momentum of about 3 GeV/c. With further increase in 
transverse momentum, our model predicts that the an- 
tiproton to pion ratio should decrease and approach the 
small value given by the pQCD. We have also calculated 
the antikaon to pion ratio as a function of transverse 
momentum. The result is similar to that for antipro- 
ton to pion ratio but with a smaller magnitude and a 
larger value at high transverse momenta. We further find 
that these ratios are reduced if antiproton production 
from coalescence of minijet partons with partons in the 
quark-gluon plasma is neglected. This effect is, however, 
comparable to that due to collective flow of quark-gluon 
plasma. To confirm the mechanism for antiproton pro- 
duction from coalescence of minijet partons with partons 
in the quark-gluon plasma thus requires a quantitative 
understanding of collective flow effects in the quark-gluon 
plasma. 

We have also studied elliptic flows of hadrons using 
quark and antiquark elliptic flows fitted to the measured 
pion and kaon elliptic flows. Predicted proton and A 
elliptic flows agree with those measured in experiments. 
We have further predicted the f2 elliptic flow, which is 
smaller than other hadrons due to smaller strange quark 
elliptic flow than that of light quarks. 

The parton coalescence model based on the test parti- 
cle Monte-Carlo method can be further extended to in- 
clude collision dynamics of partons and hadrons using 
parton and hadron transport models. This would make 
it possible to study effects due to final-state hadronic 
scattering on the transverse momentum spectra and el- 
liptic flows of hadrons, which have been neglected in 
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present study as we have compared hadron spectra el- 
liptic flows from parton coalescence directly with exper- 
imental data. Also, including parton scatterings would 
allow us to treat properly the parton spectrum around 
the cutoff momentum po, resulting in a smooth spectrum 
at this momentum. It further makes it possible to deter- 
mine both transverse and elliptic flows of partons from 
the collision dynamics instead of treating them as input 
as in the present study. Moreover, including expansion 
dynamics of the system will ensure that the total en- 
tropy does not decrease even though the entropy density 
is reduced when hadrons are formed from coalescence of 
partons. 
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